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Abstract 

In [8], Chenciner and Jimenez-Perez showed that the range of the 
spectra of the angular momenta of all the rigid motions of a fixed 
central configuration in a general Euclidean space form a convex poly¬ 
tope. In this note we explain how this result follows from a general 
“real” convexity theorem of O’Shea and Sjamaar in symplectic geom¬ 
etry [53]. Finally, we provide a representation-theoretic description of 
the pushforward of the normalized measure under the real moment 
map for Riemannian symmetric pairs. 


1 Introduction 

An n-body configuration x — (sq, • ■ • , x n ) in a Euclidean space E with masses 
mi, ■ ■ ■ , m n > 0 moving in a Newtonian force field F = 'VU(x) with reduced 
center of mass 'Yh m k x k/ 'Yh m k = 0 is called balanced with factor A if 

\ / U(x ) = —A xm 

for A : E —> E a symmetric linear operator on E and m = diag(mi, • • • , m n ) 
the mass matrix [1]. This is an algebraic equation with presumably an abun¬ 
dance of solutions for large n. It is clear that for p > 0 and k G SO (E) the 
similar configuration fikx is again balanced with factor A//i 3 . 

Let x be such a balanced configuration in E with factor A. If Z : E —y E 
is a skew symmetric linear operator and satisfies Z 2 = —A, then the rigid 
motion t (->■ z(t) = exp (tZ)x is a solution of Newton’s equation 

zm = VU(z). 
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Chenciner and Jimenez-Perez have shown that the range of the spectra of 
the angular momenta of all such rigid motions is a convex polytope [8J, which 
is subsequently used by Chenciner in the analysis of bifurcation of relative 
equilibrium motions of the n-body problem [6], 

In this note, we will show that this result is just an immediate consequence 
of a convexity theorem of O’Shea and Sjamaar in real moment map geometry 
[33], which will be reviewed in particular in the setting of a pair of real 
reductive Lie algebras. We have made an effort to to write a pedestrian 
exposition. For this reason, we have restricted ourselves to the case of central 
configurations, for which A = A is just a scalar operator. Indeed, the analysis 
of the spectra range of the angular momentum of a balanced configuration 
breaks down to this case, as has been explained in [5]. 

Finally, one may naturally ask about the density of complex structures 
corresponding to the same spectrum in the range of the spectra of the angu¬ 
lar momenta of rigid motions. An explicit description of this density requires 
more involved work and will be a question for future research. Nevertheless, 
motivated by this, and in line with O’Shea-Sjamaar’s study of the real mo¬ 
ment map, we shall give a description of the pushforward of the normalized 
invariant measure by the real moment map for Riemannian symmetric pairs. 
This provides a real version of the push-forward of the Liouvillc measure by 
the moment map [T2] in this special case. 

2 The n-body problem in Euclidean space of 
arbitrary finite dimension 

The Newtonian n-body problem in a finite dimensional Euclidean space E 
with inner product (•, •) is the study of the dynamics of n point particles 
with positions x k G E and masses rn k > 0, with time evolution according to 
Newton’s laws of motion 


m k x k = y^mjm k {xj - x k )/\xj - x k \ 3 

for k — 1, • • • , n. A configuration x = (aq, • • • , x n ) in E n is a row vector with 
entries vectors in E. Its dual configuration x* on E n is then a column vector 
with the corresponding dual vectors on E as entries. Here any vector in E 
gives rise to a dual vector on E by taking the inner product with that vector. 
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For example, with this notation x*x is the n x n Gram matrix of the position 
configuration, while xmx* is the symmetric linear operator on E sending v 
to ^2 k m k (x k ,v)x k . Here m = diag(mi, • • • ,m n ) is the mass matrix. 

The negative of the potential energy (which is also called the force func¬ 
tion by Lagrange) 

Ujx) = m.jmk/\xj - x k \ 

j<k 

is a solution of the equations 

'VkU(x) = y^/rrifmkixj - x k )/\xj - x k \ 3 

with Vfc the gradient with respect to the vector x k G E. If we denote 
Vf/(x) = (Vif/(x), • • • , W n U(x )) G E n , then the equations of motion can be 
written in the form 

x = y, ym = VU ( x ) 

as a first order system. We denote K(y) = tv(y*ym)/2 for the kinetic energy. 
The total energy is thus defined by H(x,y) = K ( y ) — U(x), and is a conserved 
quantity: Indeed, we have H = tr (y*ym) — tr(i;*Vf/(a;)) = 0. 

The total linear momentum p — m kVk G E is also conserved, which in 
turn implies that the center of mass c = s f2 i m k x k / Y^m k G E has uniform 
rectilinear motion. By the center of mass reduction we may assume that 
c = p = 0, which will be done from now on. 

For the position-velocity pair (x, y) G E n x E n , the total angular momen¬ 
tum is defined by 

L = ymx* — xmy *, 

which is a skew symmetric linear operator on E. Since 

L = ymx* — xmy* = (VU(x))x* — x(VU(x))* 


is the linear operator on E sending v to 


y: mjm k (xk, V ) 


Xj Xk 
Xj - x k \ 


V- (Xj,v) - (x k ,v) 

, mjm k — -— x k = 0 


z_/ 

j¥=k 


Xj - Xk I 


we conclude that L is conserved. The conservation of total linear momentum 
and of total angular momentum is a consequence of the Euclidean motion 
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group of E being symmetry group of the equations of motion, in accordance 
with the Noether theorem. 

For n > 3, the system is non-integrable in the sense that there are no 
other independent integrals of motion than the above, a result for algebraic 
integrals due to Bruns in 1887 [4] (substantially completed and generalized 
in [24]) and for analytic integrals due to Poincare in 1890 [37]. This work by 
Poincare on the (restricted) 3-body problem reveals the great complexity of 
the general motion in case n > 3 [ 35] . 

For n — 2 the relative position z = x± — X2 G E is a solution of the Kepler 
problem 

/jl'z = —kz/\z\ 3 z — —Xz/\z\ 3 

with k = miffl 2 , A = mi +m 2 , n — k/X > 0. For H = y\z\ 2 /2 — k/\z\ < 0 the 
motion is bounded inside the region \z\ < —k/H, and is either collinear or 
the point z moves in the Euclidean plane P spanned by z and v — z along 
an ellipse with a focus at the origin, according to the area law. Let i be 
a complex structure on P compatible with the Euclidean structure, which 
means that % : P —> P is a skew symmetric linear operator with i 2 — — 1. In 
polar coordinates (r, 6), the complex variable 

z = re ie 

is a solution of the Kepler problem if and only if (r, 6) is a solution of 

r — r9 2 = — A/r 2 , r 'O + 2 rd = 0 . 

For 6 = 0 we get the one dimensional Kepler problem r = —A/r 2 , which 
corresponds to homothetic motion of 2 in E. For r = 0 we find 9 2 = A/r 3 , 
which corresponds to rigid uniform circular motion with angular velocity 
c o = a/A/ r 3 . 

For special initial configurations x G E n , there exists initial configurations 
of velocities y G E n such that the above-mentioned Kepler orbits can be 
lifted to exact solutions of the n-body problem in E. These are the so 
called central configurations and give rise to homographic motions. They 
generalize the collinear 3-body configurations of Euler from 1767 pa and 
the planar equitriangular 3-body configurations of Lagrange from 1772 [26] . 
Later examples were found for n — 4 by Lehmann-Filkes in 1891 [27], and 
Moulton in [32], and the abundance of planar central configurations for large 
n was indicated by Dziobek, who also started to use the term “central figure 
for such a configuration” in 1899 ra- 
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Planar and spatial central configurations became a renown subject in ce¬ 
lestial mechanics, notably after the standard text book of Wintner from 1941 
[39] and a crucial paper by Smale from 1970 [38]. The question of linear sta¬ 
bility for some planar central configurations was undertaken by Moeckel in 
the eighties and nineties, generalizing the Gascheau stability condition from 
1843 for the Lagrange equilateral triangle configuration [TB]. |3(i . [25 j. [3TTJ. 
Central configurations in a Euclidean space E of arbitrary finite dimension 
were considered by Albouy and Chenciner in 1998 [TJ. We mention that it is 
not yet known to us what all central configurations are for n = 4 for arbitrary 
choice of masses, and even the finiteness problem of their number has not 
been completely settled for n = 5 (for generic choice of masses, this has been 
proven by Albouy and Kaloshin in [2]), and is yet largely open for n > 6. 
Lecture notes from 2014 by Richard Moeckel on central configurations in a 
Euclidean space of arbitrary finite dimension give a nice overview with many 
more details (also on the history of the subject), and can be found on his 
website [31] . 

3 Central configurations 

We now explain the concept of central configurations in E n and their asso¬ 
ciated homothetic, rigid and homographic motions of the Newtonian n-body 
problem in E. 

Definition 3.1. For given masses, an n-body configuration x G E n is called 
central with constant X if 

X7U(x) = —X xm 

for some scalar A 6 M. 

Since U(x ) is homogeneous of degree —1, we have 
tr(x*VU(x)) = SU(x) = -U{x) 

with £ = ^2 k (xk, ■)'VkU(x) the Euler vector held on E n , and therefore 

A = U(x)/ tr (x*xm) > 0. 

Note that central configurations are just the stationary points of the function 
U(x) under the constraint tr(x*xm)/2 = 1. Clearly, if x G E n is central with 
constant A, then for all scalars fi > 0 and all proper rigidities k G SO (E), 
the configuration jakx G E n is again central with constant A//i 3 . 
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Proposition 3.2. If x G E n is a central configuration with constant X and 
r(t) is a solution of the one dimensional Kepler problem r = —X/r 2 , then 
z(t ) = r(t)x is a homothetic motion of the n-body problem. Conversely, any 
homothetic solution z(t ) = r(t)x of the n-body problem can be expressed in 
this way for some central configuration x G E n . 

Proof. Indeed, if x G E n is a central configuration with constant A and the 
real function r{t) is a solution of r = —X/r 2 , then the motion z(t) = r{t)x 
satisfies 

zm = rxm = —Xxm/r 2 = WU(x)/r 2 = 'VU(z), 

since VU(x) is homogeneous in x of degree —2. 

Conversely, suppose z{t ) = r(t)x is a solution of the n-body problem 
for some real function r{t). By substitution into the equation of motion 
zm = VU (z), we obtain r 2 rxm = VU(x). Hence r 2 r = —A for some constant 
A G M, and so WU(x ) = —A xm, thus x is a central configuration with 
constant A. □ 

We recall that a compatible complex structure on if is a skew symmetric 
linear operator J : E —>• E with J 2 = — 1. A neccesary and sufficient 
condition for such J to exist is that E has even dimension. 

Proposition 3.3. Suppose x G E n is a central configuration in E with con¬ 
stant X = to 2 > 0. Any compatible complex structure J on E gives rise to 
a rigid motion t (->• z(t ) = exp(tuJ)x of the n-body problem. Conversely, if 
E is spanned by x, then any rigid motion solution of zm = WU(z) is of this 
form. 

Proof. Indeed, we have zm = —u 2 zm = WU(z) since z is central with scalar 
A = u> 2 . 

Conversely, if Z : E —» E is a skew symmetric operator, then the rigid 
motion z(t) = exp (tZ)x of the central configuration x with scalar A = u 2 is 
a solution of zm = 'VU(z) if and only if Z 2 x = —Xx. Since by assumption E 
is spanned by x, we arrive at Z = uJ with J a compatible complex structure 
on E. □ 

Homothetic and rigid motions of a central configuration are both special 
cases of the more general homographic motions. 
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Theorem 3.4. Suppose t (->• (r, 9) is a solution of the planar Kepler problem 

z = —A^/|^r| 3 , z = re 10 r — rd 2 = —A/r 2 , rd + 2 rO = 0 

in polar coordinates. If x G E n is a central configuration with constant A and 
J is a compatible complex structure on E, then 

t (->• z(t) = r(t ) exp (9(t)J)x 

is a homographic motion of the n-body problem. 

Proof. This is just the standard derivation of the equations of motion for the 
Kepler problem in polar coordinates. Indeed, let x G E n be a fixed central 
configuration with constant A, thus WU(x) = —A xm holds by definition. We 
have to check that 

z = r(t)e e ^ >J x 

is a solution of the equations of motion zm = VU(z) for the n-body problem. 
By differentiation, we have 

A = r _1 {(r — rd 2 ) + (rO + 2 rO)J}z, 

and by assumption, we get z = —A r~ 3 z. Since z is central with constant 
A/r 3 , we find VU(z) = —A r~ 3 zm, and hence zm = WU(z) is satisfied. □ 

Note that the term “homographic” in the terminology “homographic mo¬ 
tion” , though commonly used by celestial mechanists, should not be confused 
with the term “homography” in the geometric sense, which is synonymous 
to projective transformations. Under a homographic motion with negative 
total energy 

(r + rd) 2 / 2 — A/r < 0, 

each point particle xk G E traverses a Kepler ellipse in the plane spanned by 
{xk, Jxk} with one focus at the origin according to the area law, and all n 
point particles traverse similar ellipses. 

We end this section by showing that central configurations exist in high 
dimensions in abundance. Just take a (heavy) particle with mass M at the 
origin x 0 = 0 and a cloud of n (light) particles at positions aq, • • • , x n with 
equal masses m with Y2 Xi = 0 and with a sufficient symmetry. 
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Theorem 3.5. If G be a finite irreducible subgroup of the orthogonal group 
0(E), such that G acts transitively on the cloud • • ■ ,x n and for each 
i = 1, • • • , n the fixed point hyperplane in E of the stabilizer group Gi of Xi 
in G is equal to the line Mxj, then the configuration x = (xo, x±, ■ ■ ■ , x n ) with 
masses (M, m, ■ ■ ■ ,m) is central. 

Proof. The total force on the particle Xi is the sum of the forces expelled 
from the particles Xj for j i. Hence by symmetry this total force on Xi is 
fixed by G*, and therefore equal to — \xi for some scalar A*. By symmetry 
we have A* = Xj = A for all i,j > 1, i j- j. We can take A 0 = A as well, and 
hence we find a central configuration. □ 

An example of such a configuration is obtained by taking for the cloud the 
vertices of a regular polytope in the sense of Schlaffi ra. More generally, for 
any finite irreducible reflection group, one can take for the cloud the orbit 
of a nonzero vector on an extremal ray of a positive Weyl chamber. For 
example, in dimension 8 one can obtain such a central configuration with a 
cloud of 483840 particles with Weyl group symmetry of type Eg. But there 
are plenty of other examples, for example the minimal norm 4 vectors in the 
Leech lattice gives such a central configuration in dimension 24 with a cloud 
of 196560 particles (see [9j for explanations of these lattices). 

The planar central configurations with a regular n-gon for the cloud was 
deeply studied by Maxwell [28], and more recently by Hall and Moeckel 
113, E3- Their rigid motion is linearly stable for n > 7 in case m/M is 
sufficiently small (the larger n , the smaller m/M should be). The question 
of linear stability of these general symmetric central configurations, in case 
m/M is sufficiently small and for dimension at least 4, is completely open. 
The motivation of Maxwell for this work was to understand the stability of 
the rings of Saturn. His essay, published in 1859, was highly appreciated at 
the time, and won him the Adams prize for the year 1856. 


4 The spectra of the angular momenta 

Let x G E n be a central configuration with constant A > 0. By a suitable 
positive scaling we may assume that A = 1, which will be assumed in this 
section. For any compatible complex structure J : E E, we have discussed 
the rigid motion t i-» z(t) = exp (tJ)x of the Newtonian n-body problem. 


Note that J* = —J = J 1 , so J is both skew symmetric and orthogonal. 
The conserved angular momentum 

L := zmz* — zmz* = Jxmx* + xmx* J 

is a skew symmetric linear operator on E. The compatible complex structure 
J turns E into a finite dimensional Hilbert space (E, J ) with Hermitian 
form whose real part is the Euclidean inner product (•,•). Clearly L and J 
commute, and if we write X = xmx* for the so called inertia operator of the 
central configuration x, then 

K ■= LJ* = X + JXJ* 

is a nonnegative selfadjoint operator on (E , J). By definition, the real spec¬ 
trum of L is the spectrum of K , considered as an ordered subset of M + of 
cardinality equal to the complex dimension of the Hilbert space (E , J). 

What are the possible real spectra of L when J varies over all the possible 
compatible complex structures on El This question was posed by Chenciner, 
who conjectured it to be a convex polytope [5j, which was subsequently shown 
by an indirect argument by Chenciner and Jimenez-Perez [8] by realizing this 
real spectum range between two Horn-type convex polytopes, and observe 
that a combinatorial lemma by Fomin-Fulton-Li -Poon [15] affirms the coin¬ 
cidence of these two convex polytopes. 

The curious convexity property of this real spectrum range raised the 
question of finding for it a direct, conceptual proof, which is a question posed 
by Chenciner and Leclerc [7]. To present a direct proof of this convexity 
property, let us rephrase the question. 

Let j : E —> E be a fixed compatible complex structure. Any compatible 
complex structure on E is of the form J = k*jk for some k G 0(E), and 
therefore 

M := kKk* = (. kXk*) + j(kXk*)j* 

is a nonnegative selfadjoint operator on the fixed Hilbert space ( E,j ). Let 
us write s(E) for the space of symmetric operators on E, and write s(E,j) 
for its linear subspace of selfadjoint operators on (E,j). We consider s(E) as 
Euclidean space with respect to the trace form (Y, Z) = tr (YZ), and observe 
that 0(E) acts on s(E) by conjugation as orthogonal linear transformations. 
Note that the map 


s(E)^s(E,j), Y^(Y + jYj *)/2 
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is nothing else but the orthogonal projection of s(E) onto s(E, j). Clearly this 
map is equivariant for the conjugation action of the unitary group U (E,j). 
Therefore the question on the range of the spectra of the selfadjoint operator 
K on the Hilbert space (E, J ) as J varies over the space of all compatible 
complex structures on E boils down to the determination of the image under 
the so called real moment map 

fi:X^s(E,j), fx(Y) = (Y+jYj*)/2 

for the real Hamiltonian action of the unitary group U (E,j) on the connected 
isospectral class X = {kXk*\k G 0(E)} in s(E) of the inertia operator 
X = xmx* of the central configuration x. 

With these settings, the convexity result of Chenciner and Jimenez-Perez 
will be an immediate consequence of a convexity theorem for the real moment 
polytope of O’Shea and Sjamaar |J4] for real reductive Lie algebras. Their 
result will be explained in the next section. 


5 The convexity theorem 

The real general linear Lie algebra Qi(E ) of a Euclidean vector space E has 
the standard Cartan involution 9 : 0 l(E) —> gl(E ) given by 9(X) = —X*, 
and the corresponding Cartan decomposition 

0 l(E) = so(E) © 5 (E) 

as sum of +1 and —1 eigenspaces of 6. The commutator bracket turns so(E) 
in a Lie algebra, and s(E) in a representation space for 50(E). The trace form 
(A", Y) = tr(AW) on 01 (E) is a nondegenerate symmetric bilinear form, which 
is negative definite on 50(E) and positive definite on s(E). The conjugation 
representation of 0 (E) on s(E) is an orthogonal representation. 

Definition 5.1. A real reductive Lie algebra with Cartan involution is a 
pair ( 0 , 9) with Lie subalgebra 0 < 0 l(E) that is invariant under the standard 
Cartan involution 6 o/0l(E). By abuse of notation, the restriction of 6 to g 
is again denoted by 9, and is called the Cartan involution of g. We have a 
corresponding Cartan decomposition 

0 = 1© 5 , 1 = 0 n 50(E), 5 = 0 05 (E) 
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of q as sum of+ 1 and —1 eigenspaces of 9. The restriction of the trace form 
to g is called the trace form of q. It is a nondegenerate symmetric bilinear 
form, which is negative on t and positive on s. The connected Lie subgroup 
K < SO (E) with Lie algebra t has a representation on s by conjugation. 
Finally, we shall assume that K < SO (E) is compact, so as to exclude the 
case of quasi-periodic subgroups. The connected Lie subgroup G < GL (E) 
with Lie algebra g is a real reductive Lie group with K as a maximal compact 
subgroup. 

Example 5.2. If j : E —y E is a fixed complex structure on E then the 
complex general linear Lie algebra gt (E,j) gives, by restriction of scalars, a 
real reductive Lie algebra with Cartan involution. 

Definition 5.3. A real reductive Lie algebra (g, 6) is called complex if there 
is a complex structure j : g —* g making g into a complex Lie algebra, such 
that jd = —6j. This means that 6 is an antilinear involution of (g, j). Note 
that multiplication by j interchanges t and s. 

The complex general linear Lie algebra (gl (E,j),9) is a natural example 
of a complex reductive Lie algebra with Cartan involution. 

Definition 5.4. A real reductive Lie algebra (g, 6) with Cartan decomposition 
0 = t © s has a natural complexification (g c , 9) defined by 

0 C = 0 + *0, i = 


with Cartan decomposion 

0 c = u © p, u = t + is, p = 5 + it 

for the natural antilinear Cartan involution 9 on g c . The homogeneous spaces 
G/K and U/K are dual (in the sense of Elie Cartan) Riemannian symmet¬ 
ric spaces of noncompact and compact type respectively. Both spaces are 
different real forms of the complex symmetric space G c /K c with transversal 
intersection at the base point eK. Here G c is the complex Lie subgroup of 
GL (E c ) with Lie algebra g c (with E c the complexification of E), and with K c 
the complex Lie subgroup of G c Lie algebra t c = t + A. 

The following theorem collects the standard structure theory for real re¬ 
ductive Lie algebras with Cartan involution [22] . 
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Theorem 5.5. Let (g, 6) be a real reductive Lie algebra with Cartan decom¬ 
position g = £ © s. Any two maximal commutative linear subspaces in s are 
conjugated under K. If a < s is a fixed maximal commutative linear subspace, 
then the Weyl group W = Nk(A)/Ck(q) (normalizer modulo centralizer of 
a in K) acts by conjugation on a as a finite reflection group. Let a + denote 
the closure of a fixed connected component of the complement a° of all the 
mirrors in a, and call it the (closed) positive Weyl chamber. Then a + is a 
strict fundamental domain for the action of W on a, and likewise for the 
conjugation action of K on 5. 

Let g = t ® s be a real reductive Lie algebra with complexification g c = 
u © p as above. For X G a + we shall denote 

X = {kXk*- k G K} c s 

and call it the isospectral class of X in s. By construction, X is connected, 
and X = X D a + by the above theorem. If we denote 

X c = {uXu*;u G U} C p 

then X c has the structure of a complex manifold with real form X = X c D s. 
Moreover X c has a Kahler metric, whose imaginary part is the Kirillov- 
Kostant-Souriau symplectic form u ; on X c . The action of U on X c is Hamil¬ 
tonian with moment map the inclusion X c p. For this reason, we shall 
call the action of K on the real form X a real Hamiltonian action with real 
moment map the inclusion X H> s. 

We now have set up all the notations in order to formulate the convexity 
theorem of O’Shea and Sjamaar [33] in case of a real reductive Lie algebra. 

Theorem 5.6. Suppose (g, 9) < (g', 9') is a comparable pair of real reductive 
Lie algebras with Cartan involution. For X' C s' a fixed isospectral class the 
orthogonal projection n : X' —>■ s is clearly equivariant for the conjugation 
action of K, and is called the real moment map for the real Hamiltonian 
action of K on X'. Under all these assumptions, the intersection 

/i(X') O a+ 

is a convex polytope, called the moment polytope of the real Hamiltonian 
action of K on X'. 
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This theorem has a long history, and we shall mention just a few selected 
references. In case ( 0 , 6 *) < (g', O') are both complex reductive Lie algebras 
with Cartan involution the theorem is due to Heckman |21j . The result was 
generalized by Guillcmin and Sternberg, who replaced the coadjoint orbit 
X' of the overgroup K' by a complex projective manifold with a Fubini- 
Study metric h with a holomorphic linearizable action of K, which leaves 
the symplectic form u ; = A/i invariant, and /i the moment map for this 
Hamiltonian action of K [18] . This result was also obtained by Mumford, 
published in the appendix of a paper by Ness [33]. This is the non-abelian 
convexity theorem in the Kahler case, which generalizes the former Abelian 
convexity theorem of Atiyah [3], and that of Guillcmin and Sternberg m- 
The proof of the general case without assuming the symplectic manifold to 
be Kahler was found by Kirwan [25]. These works were all done in the early 
eighties with many more exciting developments in moment map geometry. 

It took almost two decades before O’Shea and Sjamaar discovered the 
natural real setting of the convexity theorem, which generalizes the Abelian 
real convexity theorem of Duistermaat HU- 

Indeed, consider the commutative diagram 

s' D X! -> X' c C p' 

n fj, 

5 D fi(X') -> l4X' c ) C p 

with X' c = {uXu*;u G U'}. As before, X' c can be canonically identified with 
a coadjoint orbit of U'. Therefore it has a natural symplectic form u/, for 
which the action of U' by conjugation is Hamiltonian with moment map the 
inclusion X' c > p'. The restriction of the action from U' to U gives a moment 
map fi : X' c —> p, which is just the restriction of the orthogonal projection 
p ; —> p. 

The space X! c has a natural antisymplectic involution r, which is just the 
restriction of the antiinvolution —O' of p' = s' © it', taken +1 on s' and —1 
on it'. In turn, the fixed point locus of r on X' c is just X' = X' c fl s'. Hence 
the map /i : X' —> s is nothing but the restriction of fi : X! c —> p to the 
real form X!. This explains our use of the terms real moment map and real 
Hamiltonian action. 

If b C p is a maximal commutative subspace with f) fl s = a and p + is an 
adapted positive Weyl chamber, in the sense that p + fl a = a + , then 

M(X' C ) n a+ = (/!(£') n f)+) n a 
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is a convex polytope by the convexity theorem of Heckman. 

Theorem 15.61 is therefore a direct consequence of the following result, 
which is also due to O’Shea and Sjamaar. 

Theorem 5.7. We have /r(X') fl a+ = //(£') fl a + . 

We have restricted ourselves to the case of (co)adjoint orbits for a real 
reductive Lie algebra, which both suffices for our purpose and keeps the expo¬ 
sition as concrete as possible. In their paper, O’Shea and Sjamaar formulated 
everything in the general setting of a Hamiltonian action of a connected com¬ 
pact Lie group U on a connected symplectic manifold (M,cu). Suppose that 
the group U has an involution 9 with fixed point group K, and the space 
(M,oj) has an antisymplectic involution r with M T not empty. These two 
structures are assumed to be compatible, in the sense that 

r(ux) = 9(u)t(x) and /i(r(x)) = —9(yi(x)) 

for all u G U and x G M. Linder these conditions, O’Shea and Sjamaar 
obtained the following general result 

Theorem 5.8. We have p(M T ) 0 a + = p(M) fl a + and the right hand side 

fjb(M) n a+ = n fj+) n o 

is indeed a convex polytope by the convexity theorem of Kirwan. 

It is readily seen that this implies Theorem 15.71 


6 Pushforward of the normalized measure by 
the real moment map 

We start this section by explaining the notion of Gelfand pairs, and their 
associated harmonic analysis. 

Harmonic analysis for Gelfand pairs 

Definition 6.1. A locally compact unimodular topological group G with a 
compact subgroup K < G is called a Gelfand pair if the natural unitary 
representation of G on L 2 (G/K,dx ) decomposes in a multiplicity free way. 
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It can be shown that this definition is equivalent to the following one: 

Definition 6.2. A pair K < G is called a Gelfand pair, if for any irreducible 
unitary representation (V, (•, •)) of G, the restriction from G to K contains 
the trivial representation of K with multiplicity at most 1. 

Definition 6.3. For a Gelfand pair K < G, an irreducible unitary repre¬ 
sentation (V, (•, •)) of G is called spherical ifV K = Cv has dimension 1 for 
some v G V with (v,v) = 1. The function 

G 3 g h- <p v (g) = (gv,v) 

is called the elementary spherical function associated with the spherical rep¬ 
resentation V. 

Note that elementary spherical functions are normalized by 4>v( e ) — 1- 

Definition 6.4. Any function on G that is both left and right invariant under 
K is called a spherical function. 

Yet, a third equivalent definition for a Gelfand pair is the following: 

Definition 6.5. The pair K < G is a Gelfand pair if the Hecke algebra 
TL(G/K) of continuous spherical functions on G with compact support is 
commutative with respect to the convolution product. 

The elementary spherical functions are the simultaneous eigenfunctions 
for the commutative algebra TL{G/K ), acting as convolution integral opera¬ 
tors on the space of spherical functions. 

Finally, in case that G is a connected Lie group, there is a fourth equiv¬ 
alent definition for a Gelfand pair: 

Definition 6.6. For a connected Lie group G, the pair K < G is a Gelfand 
pair if and only if the algebra V(G/K) of linear differential operators on 
G/K , which are invariant under G, is commutative. 

Similarly, the elementary spherical functions are the simultaneous eigen¬ 
functions for the commutative algebra V(G/K), acting as invariant differen¬ 
tial operators on the space of spherical functions. 

Under all these equivalent conditions, the abstract spherical inversion 
theorem gives the existence of a unique positive measure gp on the set G/K of 
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equivalence classes of unitary irreducible spherical representations of K < G, 
called the spherical Plancherel measure , such that 


cp{x) 


L 


G/K 


( p(v) ( p v (x)dMy) 


for all <p G H(G/K), with 


<KV) 


L 


G/K 


(f>(x)<j>v{x)dx 


the so called spherical Fourier transform of f> e %{G/K). 

The case that K is the trivial subgroup of G = M + or G = {z e C x ; \z\ = 
1} gives the classical inversion formula for Fourier integrals and Fourier series 
respectively. 

Harmonic analysis for Riemannian symmetric pairs 

After a brief exposition of the harmonic analysis for general Gelfand pairs, 
we now come to certain particular cases of our interest. In the notation 
of the previous section, these are the Riemannian symmetric space G/K of 
noncompact type, its compact dual Riemannian symmetric space U/K and, 
finally, the intermediate flat tangent space s, considered as homogeneous 
space for the so called Cartan motion group s x K. 

The spherical inversion formula was made explicit in the case K < G by 
Harish-Chandra [20] with simplifications by Hclgason, Gangolli and Rosen¬ 
berg (cf. [23]). Harish-Chandra enlarged the set G/K of equivalence classes 
of spherical irreducible unitary representations of G to the set G/K of equiv¬ 
alence classes of spherical continuous irreducible representations of G on a 
Hilbert space, which are only unitary for the subgroup K. He showed that 
G/K = a c /W and derived the Harish-Chandra isomorphism 


V(G/K) = Sa?,D^ lD , 


in which SaY denotes the symmetric algebra of af . The associated elemen¬ 
tary spherical functions are given by the Harish-Chandra integral formula 
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with Iwasawa decomposition G = KAN,g = k(g)a(g)n(g), Iwasawa pro¬ 
jection A(g) = log a(g), the restricted Weyl vector p (half sum of positive 
restricted roots counting multiplicities) and the normalized Haar measure dk 
on K. 

These elementary spherical functions are solutions of the system of dif¬ 
ferential equations 

Dcfx = 7d(A) 0 a , D G V(G/K) 

with normalization <fi\(e) = 1 as before. The spherical inversion theorem now 
takes the form 

* (x)= w\L kx)Mx) Ww 

with spherical Fourier transform 

0(A) = / (t>(x)(t)x{x)dx, 

Jg/k 

the Lebesgue measure fi^ on ia and the Harish-Chandra c-function A (->• c(A), 
given as an explicit product of T-factors by the Gindikin-Karpelevic formula. 

The pair K < s x K, with the semidirect product 5 xi K acting on s via 
rotations and translations, is a Gelfand pair as well, and the group s x K is 
called the Cartan motion group of the space s. The algebra T>(s) of invariant 
linear differential operators is isomorphic to the algebra Ss^ = So)! 1 of K- 
invariant linear differential operators on 5 with constant coefficients. Its 
simultaneous eigenfunctions are the symmetrized plane waves 

MX) = [ eX kXk *U,k 

Jk 

normalized by ij)\ (0) = 1 for all A G a c and X G 5. The spherical inversion 
theorem is a direct consequence of the classical inversion theorem for the 
Euclidean Fourier transform on 5, applied for functions invariant under K. 
In a sense, we can consider this theory on the flat space s as the classical 
spectral limit of the above Harish-Chandra theory for the curved space G/K, 
by the help of the following formula 

Proposition 6.7. We have 

i>x(X) = lim 0 nA (exp(X/n)) 

n—too 

for all A G a c and X G s. 
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Proof. By Harish-Chandra’s integral formula 


0 A (expX)= / e P~p^vX-k)) dk = 

J K 


e (\-p,A{exp(Ad(k)X))) dk 


'K 


we have 


0 nA (exp(A 7 n)) = / e ^-pA{eMAd(.k)X/n))) dk 

Jk 


On s, the infinitesimal Iwasawa projection 5 —y a coincides with the orthogo¬ 
nal projection s — > o. Indeed, if X G s has infinitesimal Iwasawa decomposi¬ 
tion X = Y+H+Z , for which Y G t, H G a, Z G n, then X = H+(Z—6Z)/ 2, 
which means that H is also the orthogonal projection of X on a, as we have 
the orthogonal decomposition s = a © (5 fl (n © 0n)). We therefore deduce 
that 


lirn 0 nA (exp {X/n)) = 

n—>oo 


^(Ajlimn-joo nA(exp(Ad(k)X/n))) dk _ / £ (X,Ad(k)X) dk 


'K 


11 < 


which by definition is equal to fj\(X), for all A G a, and thus for all A G a c . □ 

The elementary spherical function (\,x) 1 —> 4>\(x) is holomorphic and 
Weyl group invariant in the spectral variable A G a c , and real analytic in the 
space variable x G G/K, or, in other words, holomorphic in the space variable 
x taken from a suitable tubular neighborhood of G/K in the complexified 
space G c /K c . It has a holomorphic extension to all of G c /K c if and only if 
(A — p) lies in the intersection L fl a + of a suitable lattice L with the positive 
chamber a + C a, given in explicit terms of the restricted root system by 
the Cartan-Helgason theorem ( [231, Ch.V, Theorem 4.1]). The corresponding 
irreducible spherical representation 14(A) (with highest weight (A — p)) is then 
finite dimensional, and unitary for the compact form U of G c . If v G V(\) K 
is a normalized spherical vector, then <f>\(u) = (uv,v) for u G U with (•, •) 
the invariant Hermitian form on 14(A). 


Pushforward of the normalized measure 

After this survey of the theory of spherical functions, we can finally explain 
the meaning of the pushforward under the real moment map p, : X' —» s 
of the normalized invariant measure on X' in the notation of Theorem 15.61 
in terms of spherical representation theory. Let A' G L’ + = [K fl a' + ) + p' 
and (I4(A 7 ), (•,•)) be the associated finite dimensional spherical irreducible 
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unitary representation of U' with normalized spherical vector v' G V{X') K '. 
Let (j)\ >(u') = (u'v',v') be the associated elementary spherical function on 
U'/K'. Its restriction to the totally geodesic submanifold U/K < U'/K' is 
given by 

<t>\’(u) = 22 rn X ’(\)<f>\(u) 

A &L+ 

with my (A) = (v\,v\) if v' = J2 x v x is the primary decomposition of v' into 
components v x for A G L + of spherical vectors for the Gclfand pair K < U. 
For A' G L' + , let 

p : X' x , = {k'X'k'*] k! G K'} -G s 

be the real moment map for the real Hamiltonian action of A', and let d$' be 
the normalized K'~ invariant measure on X\,, so that / dp' — 1. 

Theorem 6.8. Let p i-> 8(p — A) be the Dirac delta distribution on a with 
unit mass at A. The probability measure p x > on a + given by 

dp x >(p) = lim m n y(X)8(p - X/n) 

n—too LJ 
\eL+ 

describes the pushforward measure p*(dp') on s by the relation 

f /(A)^(dy')(A) = f f(X)dpy(X) 

J s J a+ 

for all continuous functions f on s, which are invariant under K. 

Proof. For n G N, A' G L' + and X G s we have 

0 n y(exp (X/n)) = 22 m nx(X)<j) X {exp(X/n)) 
ag l+ 

= 22 ( n '^)4 > n X (exp (X/n)) , 

AEL-|_/n 

which in turn implies 

'lfy(X) = [ MX)dpy( A) 

J a+ 

for all X G s. Hence the desired formula for py as the pushforward measure 
p*(dp') follows from the Euclidean inversion theorem for the flat space s and 
the Fubini theorem. □ 
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This theorem generalizes the result of [21] on the relation between the 
asymptotic behaviour of branching multiplicities and the pushforward of the 
Liouvillc measure under the moment map in case (g, 6) < (g', O') are both 
complex reductive Lie algebras with a Cartan involution. In that paper, 
the convexity theorem was derived from the above theorem together with a 
simple represent at ion-theoretic property. 

Some questions 

We end this section and the paper with some questions. 

Question 6.9. For A £ L + and X' 6 L' + , does the spherical irreducible repre¬ 
sentation V(A) of (g, 6) occur as subrepresentation of the spherical irreducible 
representation K(A') of (g', 6 1 ') if and only 7ny(X) > 0? 

Question 6.10. Is it possible to generalize the results of this section to the 
general Hamiltonian setting, in line with 0 ’Shea and Sjamaar? 

Question 6.11. Is it possible to give an explicit description of the pushfor¬ 
ward measure in the example of Chenciner? 
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